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Summary 

Work  completed  under  this  grant  on  polar  fluids  includes  an  extension 
of  perturbation  theory  to  polar  molecules  with  point  octopoles  and 
hexadecapol es  as  well  as  a  preliminary  investigation  of  molecules  with 
non-spherical  cores.  In  addition,  the  ordering  of  dipolar  molecules  near 
a  flat  electrified  wall  has  been  investigated  in  the  linearized  hypernetted 
chain  (LHNC)  and  quadratic  hypernetted  chain  (QHNC)  approximations.  This 
has  led  to  a  new  study  of  el ectrostriction  and  non-linear  effects  in  the 
polarization  density  of  dipolar  fluids. 


This  report  describes  research  on  polar  fluids  carried  out  at  the 
University  of  Maine  at  Orono  under  contract  from  the  Office  of  Naval  Research 
Apart  from  the  principal  investigator  (Professor  J.  C.  Rasaiah),  a  visiting 
fellow.  Professor  Noel  Thompson  from  Australia,  a  graduate  student,  Mr.  John 
Eggebrecht  and  a  postdoctoral  fellow.  Dr.  Dennis  Isbister  from  the  Royal 
Military  College,  University  of  New  South  Wales  in  Australia,  worked  on  this 
project;  the  last  two  receiving  partial  financial  support  under  the  present 
contract. 

The  primary  objective  of  this  investigation  was  to  extend  the  theory 
of  polar  fluids  developed  by  Stell,  Rasaiah  and  Narang^  to  (a)  spherical 
molecules  with  point-octopoles  and  poi nt-hexadecapol es  which  could  serve 
as  a  reasonable  model  for  methane  and  sulfur  hexafluoride;  (b)  non-spherical 
molecules  with  dipoles  and/or  quadruples  such  as  HCN  and  0.,.  A  second 
objective  was  to  extend  our  study  of  dipoles  to  surface  adsorption  near 
an  electrified  wall  . 

A  Monte  Carlo  calculation  of  the  free  energy  changes  that  accompany 
the  switching  on  of  polar  interactions  between  the  molecules  was  planned 
as  an  unambiguous  test  of  the  accuracy  of  the  theory  when  applied  to  non- 
spherical  cores.  Some  aspects  of  this  work  are  still  in  progress,  while 
parts  (a)  and  (b)  of  the  research  described  have  been  completed.  In 
addition,  a  very  promising  and  important  theoretical  investigation  of  the 
adsorption  of  dipoles  at  a  wall  in  the  presence  of  an  electric  field  has 
been  initiated  with  the  assistance  of  Dr.  Dennis  Isbister  and  John 
Eggebrecht.  In  what  follows,  the  research  that  has  been  completed  with 
support  from  ONR  is  described  in  greater  detail  and  future  plans  to  continue 
this  work  are  outlined. 
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I.  Perturbation  Theory  for  Polar  Fluids 

While  the  potential  energies  of  interaction  amongst  point-dipoles 

and  point-quadrupol  esare  well  known,  there  are  inconsistencies  and  ambigu- 

2 

ities  in  the  literature  when  higher  order  multipoles  were  included. 

It  was  therefore  decided  to  calculate  these  energies  from  first  principles, 
limiting  our  attention  only  to  octopolar  and  hexadecapolar  interactions 
beyond  quadrupolar  charge  distributions.  It  became  necessary  therefore  to 
perform  the  tedious  contractions  of  multipolar  and  field  tensors  the 
eight  index  field  tensor  for  the  hexadecapole-hexadecapole  energy  contains 
nearly  700  terms.2  The  two-body  potentials  beyond  the  simple  dipole-dipole, 
dipole-quadrupole  and  quadrupol e-quadrupole  interactions  were  found  to  be 
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where  v,  e,  a,  and  i  are  the  dipole,  quadrupole,  octopole  and  hexadecapole 
moments  and 

cose..  =  cose.cose  .-sine. sine, cos(<t>.-4>.) 
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where  (e..,4>.)  are  the  polar  and  azimuthal  angles  of  orientation  of  the 
axis  of  symmetry  of  the  ith  molecule. 

The  potential  energy  functions  were  then  expanded  in  spherical 
harmonics,  to  allow  analytic  evaluation  of  the  angular  integrals  in  the 
x-expansion  of  the  free  energy 
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By  writing 
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where  u.  =  and  S^1  (w.)  are  the  spherical  harmonics,  the 

spatial  and  angular  variables  of  the  potential  energy  can  be  separated. 

Using  the  orthogonality  properties  of  the  spherical  harmonics,  and  the 
previously  derived  potential  energy  functions  we  have  found  the  following 
expansion  coefficients  for  the  higher  order  multipole  interactions  (t-j  =  1, 
2,  3,  and  4  correspond  to  dipole,  quadrupole,  octopole  and  hexadecapole 
moments ) 
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The  general  pattern  of  these  coefficients  allows  a  check  on  the  tedious 
and  mistake-prone  calculations. 
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With  these  results  the  terms  of  0(a  )  in  the  free  energy  expansion  can 
be  calculated.  To  go  beyond  this,  using  a  Pade  approximant  for  the  free 
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energy  in  the  manner  proposed  by  Stell,  Rasaiah  and  Narang,  requires  at 
least  the  term  of  0(x  ).  To  this  end  the  corresponding  angularly  averaged 
three-body  potentials  for  higher  multipoles  were  needed.  They  were  not 

3 

available  in  the  literature  but  were  derived  in  the  manner  of  R.  J.  Bell 
and  Rasaiah  and  Stell.4 
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Using  quantum-mechanical  perturbation  theory  Bell  has  shown  that  the 
non-additive  three-body  potential  due  to  fluctuating  dipoles  and  multipoles 
may  be  expressed  as  a  function  of  the  intermol ecul ar  distances  and  the 
interior  angles  of  the  triangle  which  they  form 


< _ i _ _ _ — - - — - - - 

Rasaiah  and  Stell  recognized  that  the  corresponding  three-body  term  in  the 
perturbation  theory  of  polar  fluids  had  the  same  form 
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where  Z£l£2£3  is  a  coefficient  determined  by  and  2.3. 

Bell  derived  the  geometrical  factor  W  (r10,r, 

dipoles,  quadrupoles  and  all  triplet  combinations  of  them, 
his  work  to  include  octoooles  and  hexadecapol es ,  we  find 
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The  geometric  terms  W  ,  W  ,  W  ,  and  W  have  yet  to  be  obtained. 
Work  on  this  is  proceeding. 

£  £  £ 

To  determine  the  specific  coefficient  Z  1  2  3  in  the  three-body 
potential  we  made  use  of  the  spherical  harmonic  expansion  coefficients  and 
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the  method  outlined  by  Rasaiah  and  Stell  to  obtain 
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From  these  sequences  one  may  predict  the  values  of  L  1  2  3  even  when  the 
geometric  terms  beyond  ^224  ^ave  not  ^et  been  °^ta'’ne<^-  We  expect  that 
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To  complete  this  portion  of  the  work  we  needed  to  evaluate  the  remaining 
four  geometric  factors  to  obtain  the  effective  angularly-averaged  three-body 
potentials  which  contribute  to  the  term  of  0(  x  )  for  all  binary  mixtures 
with  multipole  moments  through  the  hexadecapole  moment.  The  integrals  of 
the  form 

/g  i  £(  r)u-j  r]  2  doj^dri  ^ 


123^r12,r13,r23 


IU123(r12rl3r23)  dr12dr13dr23 


were  then  numerically  evaluated  using  a  method  previously  developed  by 
Rasaiah,  Larsen  and  S tell . ^  The  result  is  the  first  and  second  terms  for 
the  lambda  expansion  from  which  a  Pade  approximant  for  the  free  energy 
can  be  formed. 

In  order  to  extend  these  perturbation  methods  to  systems  with  non- 
spherical  cores  we  must  first  properly  characterize  the  reference  system. 

To  that  end,  we  proceeded  with  a  computer  simulation  of  linear  non- 
spherical  molecules.  The  objective  was  to  compute  an  ensemble  averaged 
angularly  dependent  distribution  function  in  terms  of  the  coefficients  in  a 
spherical  harmonic  expansion  of  the  distribution  function 
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From  these  coefficients  one  may  obtain,  by  the  orthogonality  of  the  spher¬ 
ical  harmonics  in  the  expansions  of  the  distribution  and  potential  func¬ 
tions,  a  very  simple  expression  for  the  term  of  0('O  for  a  linear  molecule 
with  a  quadrupol e-moment 


BXa-j 

~ir 


4irp*0^ 

5kTo5 


f  ^^22o(  r )  ”^^221  ^22  2  ^ r ^  ^ 

o 


dr 

3 


9 


where  is  the  diameter  of  one  of  the  atoms  in  the  diatomic  molecule  and  . 
is  the  reduced  density. 

Monte  Carlo  simulations  of  the  equilibrium  properties  of  a  nomonuc 1  ear 
diatomic  molecule  with  molecular  dimensions  characteristic  of  CK  were 
perfc'.iec  to  generate  an  accurate  set  of  coefficients  g  'r,  appearing 

i  ».  ill 
1  iL 

in  the  spnerical  harmonic  expansion  of  toe  radial  distr irtuier.  -unction 
g'r,  ...  Nearly  a  million  configurations  were  generated  m  calculating 
the  results  depicted  in  Fig.  1.  The  monvnts. 


wnicr  appear  in  the  term  of  0 ( - )  in  the  free  energy  decreased  stead’lv  by 

cnly2.il  over  the  last  600K  conf  igurations .  These  calculations  are 

proper lv  -ore  accurate  than  earlier  detensi nations^ of  g  (r).  even 

.  1 .  ?p. 

thougr.  tne  programs  used  were  essentially  the  same.  Tnis  is  primarily 
because  a  larger  number  of  configurations,  with  more  frecert  sar.;  linos-  at 
intervals,  were  generated.  Ir  Tar’e  I,  the  newer  estimates  o'  tne 
tori'  of  2",  • )  are  compared  with  older  date-.,  for  reduced  auadrupc.e  moments 
corresponding  the  critical  point  and  melting  point  of  Cl.,  at  a  reduced 
density  of  0.5. 

Table  I 

Comparison  of  the  term  of  0(x)  obtained  by  numerical  integration 
of  (7)  using  the  data  for  9??n(r)-  q0?,(r),  and  g?nn(r)  from  the 
independent  studies  of  Street  and  T  i  f  des  1  ey  (S&T)  <jfld  Eggebrecht 
and  Rasaiah  (E&R) 


Simpson 1 

$  rule 

Trapezoidal 

rule 

S&T 

E&R 

E&R 

S&T 

.812 

.283 

.328 

.373 

.414 

1  .4. 

.493 

.571 

.652 

.723 

i  .975 

.686 

.798 

.907 

1.006 
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A  number  of  difficulties  became  apparent  from  the  discrepancies  in  the  table. 
The  sources  of  error  in  this  calculation  were  found  to  be:  1)  the  large 
fluctuations  in  the  moments  of  g^Cr);  2)  the  sensitivity  of  the  function 
to  the  method  of  integration.  The  greater  similarities  in  our  estimates  of 
the  integral  may  be  attributed  to  the  smaller  increments  used  in  our  version 
for  the  radial  distance. 

II.  Adsorption  of  Dipoles  and  Ions  at  a  Wall  in  the  Presence  of  an  Electric 
Fi  eld 

While  our  understanding  of  the  behavior  of  liquids  has  improved 
dramatically  over  the  past  ten  years,  the  physics  and  chemistry  of  surface 
phenomena  is  still  at  a  primitive  stage.  One  of  the  outstanding  problems 
in  this  area  is  the  elucidation  of  the  density  profiles  of  ions  and  dipoles 
at  a  charged  surface;  a  detailed  knowledge  of  this  is  fundamental  to  our 
understanding  of  phenomena  ranging  from  the  behavior  of  electrodes  in 
batteries  to  selective  adsorption  and  transport  across  membranes.  It  is 
convenient  to  study  ions  and  dipoles  at  a  surface  separately  before  trying 
to  gain  insight  into  the  behavior  of  a  mixture  of  these  two.  We  have, 
therefore,  embarked  upon  a  program  in  which  the  following  problems  will  be 
studies  in  sequence. 

1)  Dipoles  at  a  wall  in  the  presence  of  an  electric  field. 

2)  Ions  at  a  charged  interface. 

3)  Dipoles  and  charges  at  a  wall  in  the  presence  of  an  electric  field. 

Investigation  of  the  first  problem  has  progressed  very  rapidly  in 

col laboration  with  Dr.  Dennis  Isbister  and  John  Eggebrecht.  Professor  Noel 
Thompson,  a  visitor  on  sabbatical  leave  from  Australia,  Dr.  Dennis  Isbister, 
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and  I  have  made  progress  on  the  second  problem,  while  work  on  mixtures  of 
ions  and  dipoles  at  a  wall  has  only  just  begun.  In  what  follows,  we  provide 
a  brief  description  of  our  work  on  dipoles  at  a  wall. 

The  effect  of  a  wall  is  obtained  by  allowing  a  mixed  system  of  at  least 
two  species  1  and  2  to  change  in  such  a  way  that  the  density  of  one  species 
(2  say)  goes  to  zero  accompanied  by  an  increased  in  diameter  R^.  It  is  the 
growth  of  this  single  particle  (2)  that  emerges  finally  as  a  wall,  and  if 
the  wall-particle  had  a  dipole-moment,  the  zero  density  -  infinite  radius 
limit  could  be  made  to  produce  an  electric  field  as  well  emerging  at  any 
chosen  angle  from  the  wall.  The  density  profile  of  each  fluid  species 
near  the  wall  is  given  by 

c,(r)  =  lim  lim  p °[h21 (p ,^2)+l ] 

where  o^cis  the  bulk  density  of  species  1,  and  h^fr)  is  total  wall-particle 
correlation  function  which  has  the  usual  invariant  expansion. 

Our  initial  investigation  of  dipoles  at  a  wall  used  the  linearized 
hypernetted  chain  (LHNC)  closure  for  the  wal 1-particle  and  particle-particle 
interactions . 

It  is  known  that  this  approximation  leads  to  fairly  accurately  bulk 
properties  for  dipoles.  We  found  that  the  contact  values  for  the  wall- 
particle  correlation  function  are  much  higher  than  those  obtained  in  the 
mean  spherical  (MS)  approximation/ 

The  linearized  hypernetted  chain  closure  is  obtained  by  taking  the 
first  term  in  the  expansion  of  the  angular  part  of  In^^i+l)  which  appears 
in  the  hypernetted  chain  (HNC)  approximation  for  the  direct  correlation 
function  .  If  the  second  term  is  retained  as  well,  a  more  accurate 


12 


description  of  the  system  results.  This  approximation,  known  as  the 

8  9 

quadratic  hypernetted  chain,  QHNC,  ’  reveals  a  new  phenomenon  that  is 
absent  in  the  MS  and  LHNC  approximations.  This  is  the  change  in  the 
density  of  the  fluid  in  an  open  system  when  the  electric  field  is  turned 
on,  and  is  called  el ectrostriction .  The  leading  term  in  the  relative 
change  in  density  at  an  infinite  distance  from  the  wall  was  found,  in 
the  QHNC  approximation,  to  be  given  by  a  term  of  0(E  ) 


k  (2)  _  Ap_ 
h  PjU 


u-i) 


2  E 


2 


Q 


where  E  is  the  electric  field,  Q  the  inverse  compressibility,  c  the  di- 

4-  2  o 

electric  constant,  and  y  =  g-  m^  6,  with  p  =  1/kT.  Approximations 
beyond  the  QHNC  theory  systematically  generate  term  of  higher  order,  so 
that  this  approximation  must  contain  the  complete  electrostriction  term 
of  0(E  )  in  the  HNC  approximation.  However,  a  discrepancy  in  was 
found  on  comparing  our  result  with  the  thermodynami cs  of  electrostriction 
in  an  open  system.^  Since  the  HNC  approximation  and  all  of  the  theories 
derived  from  it  ignore  bridge  diagrams,  the  discrepancy  in  the  term  of 
0(E  )  between  the  QHNC  result  and  thermodynamics  must  come  from  neglect  of 
these  diagrams.  The  lowest  order  bridge  diagram  was  then  calculated  and 
found  to  contribute  significantly  to  electrostriction.  Adding  this  to  the 
QHNC  result  of  0(E  ),  we  obtained  (resolving  the  above  discrepancy) 


K  _r  b(c-1)2  5(t-l)Zeyl  E2 
h  ^24-np^uy  128:rp  o  Jq 

which  shows  that  the  bridge  diagrams  cannot  be  ignored  in  calculations  of 
the  local  density  at  an  infinite  distance  away  from  the  wall.  It  appears 
likely  that  they  affect  the  density  profiles  also  when  the  dipoles  are 
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closer  to  the  wall.  This  is  an  important  result,  because  all  of  the 
theoretical  studies  of  the  ordering  of  dipoles  and  ions  at  a  wall  that 
have  been  carried  out  so  far  ignore  these  bridge  diagrams. 

Our  study  of  the  LHNC  approximation  for  the  wall -particl e  distri¬ 
bution  function  leads  to  the  constitutive  relation 

?(“>)  =  (e-1)  ?(®)/4u 

between  the  polarization  denisty  ?(<*>)  and  the  electric  field  f  with  the 
same  dielectric  constant  c  as  that  obtained  in  relating  e  to  the  two- 
particle  correlation  function  of  a  dipolar  fluid  in  the  absence  of  an 
electric  field.  The  QHNC  approximation,  however,  yields  in  addition 
non-linear  terms  in  the  electric  field  for  the  polarization  density. 

We  plan  to  extend  this  work  to  mixtures  of  dipoles  in  order  to  study  the 
selective  adsorption  of  dipoles  near  an  electrified  wall. 
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Fig.l  The  coefficients  of  the  spherical  harmonic  expansion  of  the  angle  dependant 
distribution  function  after  two  million  configurations. 
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A  theory  of  elect  restriction  which  follow1*  from  studies  of  dipolar  ordering  ai  an  electrified  wall  is  discussed  in  the  qua¬ 
dratic  hypernetted  chain  approximation,  bridge  diagrams  lor  the  wall  particle  correlation  functions  contribute  significant¬ 
ly  to  electrostriction  even  to  lowest  order  in  the  electric  field.  The  form  ol  the  constitutive  relation  between  the  polariza¬ 
tion  density  and  the  field  in  strong  fields  is  discussed. 


I .  Introduction 

The  ordering  ol' dipoles  ai  a  wall  from  which  an 
electric  field  emerges  has  been  studied  in  the  mean 
spherical  (MS)  (1  ]  and  linearized  hypernetted  chain 
(LHNC)  1 2 )  approximations.  We  icport  here  t hat  elec¬ 
trostriction  appears  as  an  added  featuie  when  the  the¬ 
ory  is  carried  beyond  the  LHNC  approximation  lot 
the  wall  particle  correlation  functions,  leading  to  a 
molecular  theory  of  this  phenomenon  m  w  hich  graphi¬ 
cal  analysis  and  integral  equation  approximations  that 
are  ubiquitous  in  the  theory  ol  fluids  may  be  exploit 
ed.  Of  the  hierarchy  of  approximations  generated  by 
the  hypernetted  chain  equation  (I INC),  we  find  that 
the  quadratic  hypernetted  chain  (01  INC)  approxima¬ 
tion  is  the  first  to  predict  electrostriction.  The  leading 
term  in  the  relative  change  in  density  an  infinite  dis¬ 
tance  away  from  a  flat  wall  is  of  0(/  *  1.  where  /■.'  is  the 
magnitude  of  the  local  field.  Approximations  beyond 
the  QHNC  theory  systematically  generate  terms  of 
higher  order,  so  that  the  QHNC  approximation  also 
contains  within  it  the  complete  electrostriction  term 
of  in  the  HNC  approximation. 


All  of  the  theories  of  electrostriction  which  have 
t hoi r  genesis  m  the  lINt  .  pproximation  ignore  the 
bridge  diagrams,  hut  we  find,  by  comparing  the  QHNC 
theory  with  the  thermodynamics  of  electrostriction 
I’j.  that  these  diagrams  must  contribute  significantly 
even  to  lowest  order  in  /  This  leads  us  to  believe  that 
they  play  an  equally  important  role  in  determining  the 
local  density  when  the  dipoles  are  closer  to  the  wall. 

#•*■  icpoit  liete  the  baicst  outline  of  our  molecular 
tlieoiy  tot  an  open  system  |4]  in  which  contributions 
Mom  the  OIINC  approximation  to  0(E*)  and  from 
the  budge  diagrams  to  the  same  order  in  the  electric 
Meld  and  tv'  lowest  oidci  in  the  density  and  the  dipole 
moment  aie  evaluated  analytically.  Our  theory  can  al¬ 
so  be  extended  to  higher  order  m  /:  but  this  will  not 
be  pursued  hcic.  The  thermodynamic  discussion  of 
Kirkwood  and  Oppenheim  [3 1  however  is  limited  to 
deriving  the  electrostriction  term  of  0(/f*),  since  it  is 
based  on  a  thermodynamic  perturbation  theory  of 
first  order  in  /. 2 .  that  uses  as  reference-state  input  the 
linear  constitutive  relation  between  the  polarization 
density  /’(“')  and  the  local  electric  field  E  for  a  field- 
independent  e. 
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Pi°°)  =  (e  1 )  E (°°)/47T.  (  1 1 

Turther  details  of  our  work  and  extensions  (<•  n 
will  be  presented  elsewhere  |4|  but  we  remark  tiuu 
oui  studs  yields  the  constitutive  relation  with  the 
same  value  of  e  as  that  obtained  from  from  relating  < 
to  the  two-particle  orientational  correlation  function 
in  the  absence  of  a  field,  and  also  provides  a  basis  toi 
determining  non-linear  terms  in  the  electric  field  loi 
the  polarization  density. 

2.  Elect  restriction  in  an  open  system 

We  employ  the  technique  discussed  in  refs.  1 1 .2 1  to 
generate  an  electric  field  by  taking  a  mixture  of  dipo¬ 
lar  hard  spheres  to  the  zero-density  limit  of  species  2 
with  an  attendant  increase  in  the  radius  AN  to  infinity  . 
under  the  constraint  that  its  dipole  moment  mi  divid¬ 
ed  by  the  cube  of  the  radius  of  the  excluded  volume 
A’,,  =A-,  +  A,  is  a  constant  A'p.  The  electric  field  / 
emerging  from  the  wall  is  related  to  A(,  by  |1 .2) 

A =  /(l(3  cos* 02  +  1 1*  l2t 

wheieci  is  a  unit  vector  dependent  upon  0 -,  jl.2|. 
the  angle  which  the  wall  dipole  makes  wirh  the  wall 
normal  n  Note  that  /■'■>  is  constant  when  the  wall  di 
pole  orientation  02  1s  fixed.  The  relationship  between 
As  and  the  Maxwell  field  E  in  the  fluid  is  of  the  tonn 

E -  |3  (2e  +  1 1)  AN  +  terms  non-linear  in  /o.  Ci 

I  lei  e  the  non-linear  terms  aie  of  magnitude  const  /  2 
+  .  when  Oj  =0,and  e  is  the  dielectric  constant  ol 

the  Hind  at  zero  field.  In  the  notation  ot  ref.  |2  | .  the 
wall  particle  correlation  function  has  the  mvanani 
expansion 

h  -i|  ( - .  A  2 . ft  ,  1  =  h 2j  (c )  +  h  let  //( 2 .  I ) 

+  /;•?,(.-)  A(2.1)  +. ...  (4i 

where  A (2, 1 )  =  *V *1  a,,d  .  1 1  =  s2' U  t  ■  s, 

in  which  .<->  and 5 1  arc  unit  vectors  m  the  directions  ot 
the  wall  dipole  and  fluid  dipole  respectively. r  is  the 
distance  of  a  fluid  dipole  (of  orientation  ft,  1  from  the 
wall,  and  U  is  the  unit  tensor.  In  this  wall  limit .  A-. 

•►oo_ 


where  A^U)  is  short -ranged  and  A  s,  is  a  constant  re¬ 
lated  to  the  electric  field  / \  (see  section  7).  Hence 

lim  li{ j  (cl  =  3A  s,  (M 

The  oiigin  ot  electiostrictioii  in  our  theo¬ 
ry  is  the  coupling  that  can  exist  between  As}  (c)  and 
Ax,  (r):  this  feedback  begins  w it  1  >  the  Ql  INC  approx¬ 
imation 

We  define  the  electros!  net  ion  ettect  A;,  as  the  rela¬ 
tive  change  in  density  ot  the  bulk  thud  when  the  dec- 
trie  field  is  switched  oil.  so  that 

K/,  =■  Ap/p}1  =  As||s*\  AN  ),  ( 2) 

where 

lAxl:.E2)  =  il  ‘J/r^t.-.As.n,  » clfl,  (M 

il  =  4rr  for  dipolar  hard  spheres  and  pi',1  is  the  density 
of  the  bulk  fluid  when  /  -  -  0.  It  A  is  the  correspond¬ 
ing  asymptotic  limit  ot'  the  angularly  averaged  direct 
cort elation  t unction  <;,  i :.  A.  il(  I.  we  find  |4|  that . 
m  the  absence  of  moleculai  poiat  iz.i,,iln  v 

A/,  =  A,  (J.  <“) 

where  <J  is  the  tnvetse  contpiessibiltty  ot  the  thud. 

(J  =  I  [pyil2\  I  , Htf. ft|  n  :  I  dll,  sin.  dz  I  101 

and  o,  j  (r.  SI  .  SI, )  is  the  driest  cone'ation  function 
of  the  bulk  fluid.  The  tliermody  uaim.  theory  ol  elec- 
trust r ict ion  |3|  tor  an  open  system  gives  A/,  toOr/'las 

A/;'  =  ( fj'Srr )  <  <»c  dp',' t  /  :  <J  (11) 

where  P  =  1  /A  A.  k  is  Boltzmann's  consiant .  7  is  the 
absolute  temperature  and  the  siipei script  (2)  means 
the  term  of  Ol/  * ). 

3.  The  OHNC  approximation  for  electrostriction 

By  considering  rite  asymptotic  limit  of o^} (c)  in  the 
invariant  expansion  tor  the  wall  particle  direct  corre¬ 
lation  (unction  c-ji  (_-.  A-,  if,  l  we  find  |4|  that  AN,  is 
related  to  /  s  in  the  (,)IIN(  apptoximat ton  by 


P  t  -i  -  itO  i  - 
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0»ixl-2: (3  cos 202  +  l)1'2  =  A21|2<?+(2A1Ip,^1  1 

+  U  (  A„p?fl?,)  3A,,/(1  +A,,)1,  (12) 

where  the  £)„(  )  functions  are  defined  itt  ref.  |2|. 

and  the  constant  A  u  is  related  quite  generally  tie.  not 
just  m  the  QMNC  theory)  to  the  fluid  dipole  moment 
m  |  by  |4  6| 

3v  =  C>+(2Aup‘iA|i)  QA-  Anp?K-],).  U3I 

whcie  f  =  4rrnijp([(3/9.  The  dielectric  constant  c  of 
the  tlmd  is  also  quite  generally  given  by  |('l 

e  =  tM2A11p?«J,»  CLl-  A„p?*?,).  (Mi 

We  also  find  from  the  z  -*  00  limit  of  the  QHN(  ap¬ 
proximation  for  c'2j  (:.  f2 .  £2} )  that 

A,1."1  =4  A 3j (3  cos202  +  I)  (15) 

and  on  using  (2),  (3),  (9)  and  (12)  ( 14)  we  have 

A;1,"1  =  ( (3;24ffP] »’)(e  -  l)2  E~/Q.  (lot 

where  the  superscript  (2)  means  again  the  term  of 
CM  1  -).  Note  that  in  (16)  A/,2’  is  independent  ol  the 
inclination  0-,  of  the  wall  dipole.  Consistency  with  the 
thermodynamic  theory  (11)  demands  that 

pV  *e  V,1  =(e  U“/3 v.  (IT) 

This  relation  is  indeed  satisfied  by  the  simple  Debye 
equation  (e  -  1  )He  +  2)  =  r.  but  the  more  exact  ex¬ 
pression  [7.8]  (for  arbitrary  y  and  p^1  -»0) 

ic  1)  (e  +  2)  =  y  ,fJ  +  ...  (lb) 

introduces  a  discrepancy  of  0(>,2)in(17).  Since  higlrer 
order  approximations  derived  from  the  I1M  equation 
do  not  produce  additional  contributions  to  A;,  ol 
0(1-  ),  the  apparent  inconsistency  lies  in  ignoring  the 
bridge  diagrams. 


4.  Bridge  diagrams  of  0(£2 ) 


The  bridge  diagram  of  0(A2)  and  to  lowest  order 
in  the  fluid  density  can  be  represented  graphically , 
when  r  by 


«‘<2.lt  = 


(l‘>) 


where  c - -  -  •  =  3A2t  D(Z.  i).  and 

■>  “ 


is  the  three-particle  direct  correlation  function  r3( A" j . 
V3 ,  A’4 )  for  the  fluid  with  .V,  =  (r/,£ l(  I.  /3* (2.  1 )  con¬ 
tributes  to  the  wall  -  particle  direct  correlation  func¬ 
tion.  In  (19)  the  open  circle  2  is  the  wall  dipole  root 
point,  the  half  black  circle  »  signifies  angular  integra¬ 
tion  over  the  orientations  of  the  bulk  fluid  dipole  1 
and  the  black  circles  (field  points)  icpresent  spatial 
and  angular  integrations  of  tlmd  dipoles.  To  evaluate 
(19)  to  lowest  order  in  pj’  and  Wj  we  use  the  low-den¬ 
sity  limit  in  terms  of  Mayer  f  functions 

r3(Y, .  *3.  -V4  )  "  /(*, .  -V3  )  /(*3,.Y4  )  f(X4.Xx ) 

(20) 

and  retain  only  the  n  =  1  term  in  the  perturbation  ex¬ 
pansion  [8.9) 

/<*,.*,)  =/n(z„) 

+  [1  +/n(r„)l  2  [0m]DU.j)irh"/(\\  (21) 

n  =  1 

for  the  Maver  f  function,  where /n  (r,,  I  is  the  corre¬ 
sponding  Mayer  function  for  the  reference  hard-sphere 
system.  The  bridge  diagram  is  now  evaluated  analyti¬ 
cally  using  llankel  transforms  and  we  find  eventually 
that  when  r  -* 

j-f#(£  1 )-  p v/:2  np\'.  (22) 

in  which  the  contribution  from  the  three -particle  cor¬ 
relation  function  appears  as  the  tlmd  virial  coefficient 
for  hard  spheres!  Adding  tins  to  ( 16).  after  using  (9), 
we  have 

^(2)  J'jHe  ->2  .  U2ftv 
h  L  24jtp(1'.v  1 28np9 

where  —  emphasizes  that  our  calculation  of  the  bridge 
diagram  of  0(1'-)  is  correct  only  to  lowest  order  in 
p^  and  bi| .  This  is  sufficient  however  to  test  the  mo¬ 
lecular  theory  against  thermodynamics  using  (18)  for 
the  dielectric  constant.  Instead  of  ( 1 7)  we  now  need 

p?3e/dp?~-(e  1)2  '3v  £  (e  Il2r.  (24) 

which  is  consistent  with  1 18)  to  0(  r2).  We  have  thus 
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shown  t hat  the  bridge  diagrams  make  significant  con¬ 
tributions  to  the  density  when  and  it  appears 

likely  that  they  affect  the  density  profiles  also  when 
the  dipoles  aie  closer  to  the  wall.  Theoretical  studies 
of  the  ordering  of  dipoles  and  ions  at  a  wall  that  have 
appeared  so  lar  ignore  these  bridge  diagrams. 


5.  The  polarization  density  f'i )  in  an  open 
system 


From  the  definition  ol  Pyz.  )  in  the  grand  ensem¬ 
ble  [ts|  we  have 

P(:.E2  i  =  (w,p'/  <2)  (25) 

In  our  discussion  Fi  is  independent  of  z  [see  (2)| . 
Using  the  invariant  expansion  for  ^t  >• we 

find,  after  doing  the  necessary  angular  integrations, 
that 

P{°°.  A\ )  =  "ijp'i'A  I. i  cos*/),  +  n*  *  e,.  (2o) 

w  hich  is  independent  of  any  approximation  for  the 
wall  particle  correlation  function!  Employing  the 
QHNC  approximation  lor  K ,j  given  in  (12), 

/’(«'. f  ;>  =  |2(.)+(:Aup^«J,) 

+  i>  (  A  , ,  pVA’ii  >  -’AV<  1  +  K/,  )|  1  (271 


=  l<* 


X 


1)  47r]|5f  2  t2e  +  1)| 

Kh  e  1 ]  ' 
.'ill  +  Kh)  2e  +  1J  ' 


(28) 


When  Kj,  =  0  (i.e.  m  the  absence  of  electrostriction  as 
in  the  MS  and  l.l  INC  approximation)  and  when  all  the 
higher  coefficients  of  /:' i  (rt  >  5)  in  (3)  arc  assumed 
zero,  we  recover  the  constitutive  relation  ( 1 ).  Con¬ 
versely  when  k tl  and  the  higher  coefficients  in  (3)  arc 
not  zero,  non-linear  terms  in  the  polarization  density 
will  appear.  Substitution  of  ( 1 5)  and  (9)  in  (28)  follow¬ 
ed  by  expansion  of  the  denominator  draws  out  the 
term  of  0(/ :  2  )  in  the  QHNC  approximation  for  the  po¬ 
larization  density  which  is 

^UllNr  ,=  f<c  1)  4rr|[3£2/(2e  +  D| 


X 


le  1 
2e  t  I 


+  0(p:52). 


(29) 


l‘>2 


I)  the  non-linear  terms  in  (3)  are  neglected 
/UIINC,oc-/-  ,=  l,f  D  4n|/ 


1  ♦ 


24 rp»y- 


(e  1)? 

2f  +  1 


I- 

L> 


+  Ol  / 5 1. 


(31 1 1 


A  more  complete  calculation  of  this  term  in  l he 
QHNC  approximation  would  t ecpiiie  the  delcmuna- 
lion  ol  i lie  coefficient  of  the  /  ,  term  in  ( 5  i  when  >■  - 
~  0.  An  exact  calculation  of  the  term  ol  OU  ' )  w  ould 
reijune  in  addition  that  the  contributions  ot  the  bridge 
diagiams  lo  eleetiostriction  |see  i22l|  and  to  the  rela¬ 
tionship  betw  een  K  ^  and  /  ,  be  meluded. 
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The  adsorption  of  dipoles  at  a  wall  in  the  presence  of  an 
electric  field:  The  RLHNC  approximation50 
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The  adsorption  of  dipolar  hard  spheres  in  the  presence  of  an  external  electric  field  has  previously  been 
studied  within  the  context  of  the  mean  spherical  approximation.  In  order  to  quantify  the  significance  of  the 
physical  trends  found  above,  the  problem  is  solved  within  the  higher  order  closure  rules  afforded  by  the 
linearized  hypernetted  chain  approximation.  I  xpressions  for  the  reduced  dipole  moment  and  the 
electric  field  strength  are  derived  using  only  the  asymptotic  forms  of  the  direct  correlation  functions  It  is 
found  that  the  favorable  orientational  correlations  between  ihe  dipolar  hard  spheres  and  the  wall  are 
underestimated  by  the  mean  spherical  approximation.  This  is  emphasized  in  the  enhanced  adsorption  of  the 
dipolar  species  (at  the  wall  itself)  for  dipoles  oriented  close  to  the  direction  of  the  field.  However,  the 
nonphysical  features  of  the  mean  spherical  approximation  (manifested  in  the  negativity  of  the  density  profile) 
are  not  fully  rectified  by  the  use  of  the  linearized  hypernetted  chain  approximation. 


I.  INTRODUCTION 

The  adsorption  of  dipolar  molecules  at  a  wall,  in  the 
presence  of  an  electric  field,  is  of  interest  in  the  study 
of  electrode  and  membrane  phenomena.  Here  the  ad¬ 
sorption  phenomenon  is  delineated  by  the  distribution  of 
molecules  at  a  particular  orientation  12,  and  distance  z 
from  a  hard  planar  wall  p,(z ,  E2,  12,).  The  electric  field 
E2  emanates  from  this  wall,  the  declination  of  the  field 
with  respect  to  the  wall  being  allowed  by  its  noncondue- 
tive  properties.  Isbister  and  Freasier1  have  investi¬ 
gated  this  problem  for  hard  dipolar  spheres  against  a 
hard  wall  using  the  mean  spherical  approximation 
(MSA).  Their  results  for  the  density  profile  p,U,  Eg,  f2,) 
of  dipoles  are  of  great  interest  even  though  they  suffer 
from  the  defect  that  the  wall  particle  density  profile 
p,(z.  E2,  S2 ,)  assumes  negative  values  at  certain  relative 
orientations  of  electric  field  (E2)  and  dipole  moment  m, 
of  the  particles  in  the  fluid  (the  dipole  orientation  is  de¬ 
noted  here  by  S2,).  However,  the  argument  leading  to 
the  electric  field  at  the  wall  is  not  swayed  by  the  approx¬ 
imation  used,  and  may  be  employed  with  more  accurate 
theories  such  as,  for  example,  the  linearized  hypernetted 
chain  (LHNC)  approximation.  While  these  theories  could 
be  expected  to  produce  better  results,  they  do  suffer  from 
the  necessity  of  employing  numerical  methods  to  a 
greater  extent  than  is  needed  to  determine  the  density 
profiles  in  the  mean  spherical  approximation. 

This  paper  is  devoted  to  a  study  of  the  wall-particle 
density  profile  using  the  linearized  hypernetted  chain 
approximation2  for  the  wall  particle  and  particle-parti¬ 
cle  interactions,  except  that  the  effects  that  are  inde¬ 
pendent  of  the  orientations  of  the  electric  field  and  the 
fluid  dipole  are  treated  exactly.  In  our  study,  these 
are  the  interactions  between  the  hard  cores  in  the 
fluid,  and  also  the  interactions  between  these  cores  and 
the  hard  wall.  We  call  this  the  renormalized  linearized 


*’ Extracted  in  part  from  the  Nl.S.  (Chemistry)  thesis  of  .1 
Eggebrecht,  (  niversltv  of  Maine  (1980). 

"Present  address:  Department  of  Chemistry,  Faculty  of 
Military  Studies,  I'niversity  of  New  South  Wales,  HMC, 
Duntroon,  ACT  2600,  Australia. 


hypernetted  chain  approximation  (RLHNC)  after  the 
nomenclature  introduced  by  Stoll  ,>nd  Weis.3  Our  re¬ 
sults  for  this  theory  are  an  improvement  over  the  cor¬ 
responding  MS  approximation  when  the  system  is  char¬ 
acterized  by  (a)  a  reduced  fluid  density  p*  -p,RiiOf  0.  573. 
(b)  a  reduced  dipole  moment  m*  ■  m,  v  liTR\j  of  v  0.5 
(or  equivalently  a  reduced  temperature  T*  s  1  n/*2  of  2), 
and  (c)  a  reduced  external  electric  field  E *  =  EZR\ ,  » u, 
of  8/3.  In  contrast  to  the  MS  approximation,  the 
RLHNC  wall-particle  density  tunctions  p^z,  E2,  f2,)  are 
only  marginally  negative  near  the  wall  for  a  dipole  ori¬ 
entation  in  direct  opposition  to  the  field  (see  Fig.  2). 
These  functions,  however,  can  become  negative  over  u 
larger  distance  z  from  the  wall  when  the  reduced  dipole 
moment  is  increased  to  1.0  without  altering  the  reduced 
electric  field  or  reduced  density  (Fig.  5).  This  sug¬ 
gests  even  higher  order  terms,  beyond  the  RLHNC  ap¬ 
proximation,  must  be  included  in  the  theory  when  the 
dipole  moment  «/,  and  the  external  electric  field  Ez  are 
both  large. 

II.  GENERAL  THEORY 

The  technique  of  producing  a  wall  next  to  a  fluid  bj 
taking  the  limiting  behavior  of  a  binary  mixture  (with 
densities  p,,  p2  and  radii  R,,  R2)  detailed  by 

lim  lim  (2.1) 

^2*  •  °7.~  0 

is  well  known.4  Isbister  and  Freasier1  have  extended 
this  to  introduce  an  electric  field,  as  well,  by  consid¬ 
ering  the  corresponding  limit  for  a  dipolar  mixture  un¬ 
der  the  restriction  that  the  dipole  moment  w2  of  par¬ 
ticle  2,  which  eventually  becomes  the  wall,  divided  by 
the  cube  of  the  radius  of  the  excluded  volume  R2,  =R2 - R, 
is  a  constant  (hereafter  called  E0): 

lim  »«*/*!,  =  E0  .  (2.2) 

8  2"  “ 

In  taking  these  limits,  in  the  specific  order  pt-  0,  R2 
— 00 ,  the  volume  of  the  system  is  allowed  to  grow  faster 
than  R2,  keeping  p,  constant  through  the  constraint  that 
p2R2—  0. 

The  magnitude  and  direction  of  the  electric  field  E2 
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I'M.  I.  Ik'K!  ...  :i\  iix  wall  ;-l  >ticd  against 

:':u  itKlinai i->n  tr  .  •*  ‘Vail  -li;  ,ilc”  v. 1  u-  '  rrsi-les  at  minus 
inlim:  an«!  pro<iact*>  * ht*  oUvlru-  ticl'i  £  .  iV.t  an^le  f*,  is  the 
inclination  ot  a  tlij-oU*  m  the  ihu-.I. 


follow,  when  this  limit  is  applied  to  the  dipole-dipole 
interaction  energy 

r21<r2l.  D2,  &,)  -  -  '-^r1  D( 2.  1)  .  (2.  3) 

r  21 

Here 

/)(2,  1)  -  •  (3P21  r,j  -  U)  •  <2  ,  (2.4) 

.C,  and  s2  are  umc  rectors  in  the  directions  oi  the  dipole 
moment  vectors  nq  and  m2.  respectively.  (3r-21r2i-U) 
is  the  dipole-dipolo  interaction  tensor,  in  which  r21  is 
a  unit  vector  in  the  direction  along  the  line  joining  par¬ 
ticles  1  and  2.  and  U  is  the  unit  tensor.  The  potential 
may  be  written  m  terms  ot  the  electric  field  E2,  pro¬ 
duced  at  the  location  of  particle  1  In  the  second  parti¬ 
cle  2. 

V’21(r21.  •  E2  .  (2.5) 

w  he  re 


In  Eq.  (2.6' 

e2  (3r21  r21  -  U)  •  x2  -  (3r21  cose-  -  -,)  (2.  7) 

and  t* i  is  the  angle  which  the  di()ole  embedded  in  particle 
2  makes  with  r21  (see  Fip.  1).  Since  the  magnitude  of 
e2  is  (3  cos2b2  -  l)1  2. 

E2 (3  cos2*?,  *  l)1 '  *■  rz  ,  (2.8) 

r  21 

where  r2  is  a  unit  vector  in  the  direction  of  E2.  Changing 
variables  to  r21  -  Kz  •  r .  and  taking  the  wall  limit,  the 
electric  field 

E2  ;£0(3cos!e2- l)WJr2  .  (2.9) 

which  shows  that  E2  is  independent  of  the  distance  z  from 
the  wall,  but  that  its  magnitude  is  determined  by  the 
strength  (through  E0)  and  the  direction  (through  S2)  of 


llu-  "wall  dipole."  which  has  receded  to  a  distance  -  ' 
from  the  wall.  At  that  distance,  the  vector  r£1  become 
perpendicular  to  the  wall,  and  the  electric  held  which 
makes  an  angle  o,  with  respect  to  this  normal  is  given 
by  the  appropriate  solution  to 


cos  o  | 


2  cose. 


(3  cost,  •  1) 


(2.  101 


The  solution  to  this  equation  gives  the  direction  oi  the 
electric  field  uniquely,  in  terms  of  the  orientation  m 
tiie  wall  dipole  al  minus  infinity  (see  Fig.  1). 

The  density  profile  E2.  S2,l  of  dipolar  molecule- 

at  a  distance  .  trom  llu  wall  depends  on  I  Is  dipole  - 
entation  CIj,  n  the  direction  ot  the  electric  tieln  E;  ano 
on  oilier  parameter.-  of  the  system.  It  :s  obtained  iron, 
the  relation1 


.  E„  h.) 


lim 

Hr  -  ‘ 


lull 


S-',l  •  1 


(2.  11' 


where  , > ,  is  Uie  bulk  density  ot  species  1  and  (o,h?!.  I.,, 
f 2 , )  is  the  total  correlation  (unction  ot  sptc  ics  1  .,no  2 
a  binary  mixture.  The  latter  is  the  solution  to  me  Urn- 
stem-  Zernike  relation 


,  1  V  , 

'21-  (*21  *  >■*,  ('? 

S  V  3  1 


(2.  12) 


where  s  S  r/r2r/D3  is  a  convolution  involving  spatial  and 
angular  integrations.  Assuming  the  invariant  expansions 


hnb 2i.  S(2.  «,)  -/(!,(>■*,)  -  Ji£,<»"„)  M2,  1)  -  //?,(.-„)/>( 2.  1)  . 

(2.  13) 

c2i(r2i>  h2.  h  j)  -  r2j(r21)  *  ("2i(r2j)  A(2.  1  -  -  c’2](r21)/i(2.  1'  . 

12. 14) 

and  using  Werthoim’s  multiplication  table"  in  Fourier 
space  for  the  angular  integrations,  the  Ornstcin-.  crabs 
relation  reduces  to  three  equations 


iS, 


V> 

'  21 


>V 

3  £ 


,  b'K' 2, 


-  i’i  ■ 


(2.  H" 


lit 


,>,(2(4 


1  ,  i 


(2.  17) 


when  *  !  dtJ  is  a  convolution  involving  only  spatial  in¬ 

tegrations.  in  Kqs.  (2.  16)  and  (2.  17) 


i£e<»")  ) -3  f  ,/s  s-' l&ts)  12.18) 

and  its  inverse  is 

fT  i'‘‘e(s)  s‘,/s  ,  12.19) 

,  *  u 

with  similar  expressions  for  r£a(r)  and  c£e(>)  in  which 
a,  d  =  1  or  2.  Inside  the  hard  cores 


-  1  , 

Kab')  =  htfb'l  -  0  ,  r  •  fia9  =  fi0  /?„ 
so  that  from  Eq.  (2.  18) 

l'a*(r)  -3 A’a,  .  f  Hae  , 

where 


(2.  20) 


(2. 21) 
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Kat~-  I  h°aa( s)s'l<ls  .  (2.22) 

•SjD 

Following  We rthei ni. 3  the  relations  for  /;£,(>)  and 
/i2l(r)  may  be  uncoupled  by  taking  linear  combinations 

'»;,(»•)  -  [/<“„(»-)  *  '■  ^fl(r)  I  3K0„  .  (2.  23) 

-  !C(r)  -  *£„(,->  |  3A-ofl  .  (2.  24) 

with  analogous  expressions  for  r*a(r)  and  c"e(r).  The 
closure  conditions  (2.20)  when  applied  to  Eqs.  (2.231 
and  (2.24)  are  equivalent  to 

li‘ai(r)  -  -  1  .  r-/f0,  .  (2.25) 

On  using  the  linear  combinations  l:‘asi(.r)  in  Eqs.  (2.  16) 
and  (2.  17),  one  finds 

**,<*•>  -  r|,(r)  -  £  A„  it  .  <•;,  .  (2.  20 

1=1 

where  p, *  -  2p,,  p"  — ,g.  Equations  (2.  15)  and  (2.26'  re¬ 
semble  the  Grnstein-Zermke  relations  for  a  binary 
fluid,  in  which  the  molecular  interactions  are  spherical¬ 
ly  symmetrical. 

In  the  limit  ft,-  -i ,  Eq.  (2.  19)  together  with  Eqs. 
(2.23)  and  (2.24)  yields 

/i?,U)  =  [2/i;,U).//ilU)*3|A'!1  .  (2.27) 

As  p2- 0,  Eqs.  (2.  15)  and  (2.26)  reduce  to 

/tfi(r)  =Cj](r)  -  pj  hi,  *  <‘i j  (2.28) 


h\\(r)  -  f21(r)  »  A' n  p*,  ft 2i  *  <•*„  . 
which  can  be  written  in  bipolar  coordinates  as 


mined  by  the  interactions  between  the  particles  in  the 
bulk  fluid.  An  analogous  equation  can  be  written  for  the 
angularly  independent  part  of  the  wall-particle  total  cor¬ 
relation  function 

h'a{:)  -  <•;,(.-)  •  2up,|z  [«-V)  -B‘U)|-  [Z35(- )  -V'U)\ 

-  P,yV>i|1(  v)lft5K)-Bs(U-v|)||  ,  (2.35) 

where  li'(; '  and  D‘{z)  have  definitions  which  correspond 
exactly  to  li‘(z)  and  D‘(z). 

In  Eqs.  (2.34)  and  (2.35),  1  cS,U),  <•?,(*)'  and  ,-’zA:  . 

<  -  j  t  ( - ) ,  are  sets  of  two  different,  but  consistent,  direct 
correlation  luiictions  with  corresponding  closures  (dis¬ 
cussed  m  the  next  section)  for  the  wall-particle  and 
particle-particle  interactions,  respectively.  The  func¬ 
tions  C|,(r)  and  rj,(r)  for  the  bulk  fluid  particles  are  ob¬ 
tained  in  an  independent  calculation  from 

A’ii<j  )  -  <  i',(r)  - /»i /rfi  *  cJi  .  (2.3t.i 

(Tn(r)  -  c'*i,(r)  ->  ft  up.  Ii*u  *  <’u  .  (2.2< 

winch  an  the  one-component  analogs  ol  Eqs.  (2.28)  ..no 
(2.  29)  first  derived  in  a  seminal  paper  by  Werlhein;. 
These  equations  carry  their  own  closures  for  /in(r)(>  /•: ., . 

and  <•,,()•)()•  Kn).  We  do  1101  actually  solve  Eos.  (2.  36' 
and  (2.35)  since  nearly  exact  results  are  available  from 
the  work  of  Verlet  and  Weis6  for  hard  spheres  and  from 
the  study  of  Waisman.  Henderson,  and  Lebowitz1  for 
hard  spheres  against  a  wall.  Our  RLHNC  approximation 
(See.  IIP  implies  that  Ii5,(c)  is  the  exact  wall-particle 
total  correlation  function  for  hard  spheres  against  a 
hard  wall.  Equation  (2.37)  has  been  solved  by  Wei  - 
iheiui.  in  the  mean  spherical  approximation. 6  w!ui< 
Patey  and  his  colleagues  have  treated  it  in  the  l.!!\C 


nr-  t  c  *  .r*(  *  »ut.  i  me  nan  w  vaivu  *t  hi 

h\ i(»')  -  *’«(»')  ~ — liiil  j  i ft  ihl^i)  j  i/ssr'„(.d  .  and  QHNC  approximations. 6  using  the  Verlet -Wei 

1  ‘n  lr~f!  ,u.  i...  ..a  ...  tl'c  , tl.r. 


(2.30) 

On  substituting  r  -  ft,  •  ; .  I  -  ft,  -  v  and  on  taking  the  wall 
limit/?,-*,  we  find,  for  ;  '0. 

h'nU) -cltU) j  ,/viV21(v)  /  r/s  s  <•*„(.«)  . 

-I i-vl 

(2.31) 

where  i  is  the  distance  of  the  center  of  the  dipolar  hard 
sphere  from  the  wall.  The  integral  between  the  limits 
-  *  and  *  may  be  simplified  by  noting  that  /i2,(  v)  -  -  1 . 
when  -«<y<  0,  and  in  addition  u-  -  yl  -i  -  v.  when 
z  >  0  and  - 00  <  y  ■-  0.  On  changing  the  order  of  integra¬ 
tion  between  -*>  and  0. 

I  ilx  >i\\(  v)  f  ils  s  ('5,(n1  I  <ls(z  -  s)  sc'n(s)  . 
a—  -'u-vl 

(2.32) 

Defining  the  functions 

B‘U)=  f  rl,(s)s,/s  .  !>'(:),  I  *  r*„(s)  s*rfs  .  (2.33) 

*0  -'ll 

the  wall-particle  ±  equations  become 
l'hU)=c\tU)+  2tr  *,,*>*,{,-  !ft*(»  )-/?*(--)  | -/>'(:)  ■■ 

•2  f  ./V  //*,,(  \  I  |ft*(*)  -  B‘(\z  -  v  )|  . 

-'ll 

(2.341 

where  the  functions  /?*(.  )  and  /)’(/)  are  entirely  deter¬ 


orv  for  tin  hard-sphere  interactions.  We  are  then  ton 
left  with  t ho  necessity  of  solving  only  Eq.  (2.34'  ur.dt  r 
the  appropriate  closures  for  the  wall-particle  and  par¬ 
ticle-particle  direction  correlation  functions. 

The  constants  Kn  and  A'21  determine  the  dipole  mo¬ 
ment  w  j  and  the  electric  field  E2  through  the  relatic>",s 

TjyV-1  -  (AlSA'nPifth)  -  V.(-  A'uPift'u'  (2.3s1 


^ y-  "  A  2|  |  2(3.(2AI1p,ft  u)  s  </.(  —  A„ii,  R  ■ , '  ,  (2.39 

where  /.  and  T  are  the  Boltzmann  constant  and  absolute 
temperature,  respectively,  and  <V..(Pi  ftj ,)  is  defined 
by: 

<i'..ti>i ftfi'  1  -  4'pt  I  cf,(r,  p,)  »•*</»•  .  (2.401 

-'o 

Equation  (2.38)  has  been  derived  by  Werthein  in  the 
mean  spherical  approximation,  but  its  extension  to  im 
l.HNC  tui  HEHNC)  approximation  is  straightforward: 
nevertheless  we  present  it  for  completeness  and  as  a 
prelude  to  the  derivation  of  (2.39).  which  Isbister  and 
1  rt  asier1  discussed  in  the  mean  spherical  approxin..: 
turn.  The  derivation  of  Eq.  (2.38)  rests  on  the  asvr:.;i- 
toll,  forte  of  c|>1(j)  (see  Eqs.  (3.1)  and  (3.6' 
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«•?,('•>  -  li{2.  1)  .  .1.-  r-  •  (2.41) 

On  inserting  tins  m 

.*?,(.')  -.■?,{,•) -3  /  ./>  .  (2.42) 

one  sees  that  rfjO)  is  short  ramanl  ami  tends  to  zero  as 
r  —  1 

as  )  -  .  (2.43) 

When  the  inverse  relation  n*f.  Kqs.  (2.  19)  and  ( 2 .  1 8)  | 


■a  hie!:  \>  uimtir.i!  It)  Kq.  (2.39).  11  should  he  t  inph.i- 

si/ed  that  I  q  .  t  2.  3b)  and  (2.39)  are  relations  dene* 
lrnjj]  tile  a.-  >  inptot  ic  tonus  of  rjj  and  e[,  m  the  Ms  .*:i  : 
I.HN'C  approximations.  The  dipole  moment  >r. ;  and  t:.« 
ideetrie  lield  K2  are  determined  Irom  An  and  A‘2I,  an* 
the  solution*  t< j  Kq.  (2.37)  have  been  obtained. 

Tile  ennsl.int  A'21  also  retrieves  the  eoetlieient.s  /  d, 
and  //21(.  ),  -A inch  appear  in  the  invariant  expansion  o \ 
the  wall -particle  correlation  iunctions 

e;,  i,y. )  •  //?,(.-) m2.  i)  -  /.•j,(.)  M2.  i)  (2. : 


-  4  |  (2.44) 

)  Jr 

is  taken  to  the  limit  >  -  x  ,  one  ..Iso  tuuis  that 


•Jh 

kT 


■3  f 


(  j i(s )  s*tis 


--3 A’ , ,  j  [2rJ,(s.  2/\,j,i.)  •  i'7,(.s,  -A'uP,)! s*ds  . 

(2.45) 

where  we  have  used  the  analogs  of  Kqs.  (2.23)  and 
(2.24)  applied  to  r?|(s)  in  the  last  stop.  The  result 
given  m  Eq.  (2.38)  follows  immediately  on  applying 
Eq.  (2.40). 


The  relation  between  A'.,  and  £  ,  is  also  derived  from 
the  asymptotic  from  ot  cf,(r)  in  the  MS  and  LHNC  ap¬ 
proximations  [see  Eqs.  (3.1)  and  '3.6)  in  Sec.  Ill |.  The 
equation  for  dipolar  mixtures  which  corresponds  to  Eq. 
(2.45)  is9 


lim 

r2l-- 


»‘i»!  I  3A'2i  ( 

7x77  =  -  )»n>  / 

r21'“  >  21  So 


l2<-2l(')s  fi,(sl|ssrfs  . 


(2.46) 

where  we  have  not  canceled  the  r21  in  the  denominator, 
because  we  intend  to  take  the  wall  limit.  Substituting 
i-ji  = /?2  - and  taking  the  limit  Rs  - *,  reduces  the  left 
hand  side  of  Eq.  (2.46)  to  n/jET,,  >~T.  and  replaces  the 
upper  limit  of  the  integral  in  Eq.  (2.46)  by  * .  The 
Fourier  transform  of  Eq.  (2.29),  which  is  an  Ornstoin- 
Zernike  equation  for  mixtures,  as  p:-0,  yields 


?»(*)  =  [i-a'„p1  ?;,(*)  |/1|, ao 


(2.47) 


from  which,  when  k-0,  we  have 


f  c‘z](s)szds  = 
"/0 


[i-A'llP?  ?;,((» 1 1 


s2rfs  . 


(2.48) 

In  the  above,  the  "'"represents  the  three  dimensional 
Fourier  transform 


t(k)  =  Jdrexpiik- r) /( i  r  | ) 


and  t*c  or  h. 

Since  /i5,(s)  is  a  short-ranged  function  equal  to  -  1  for 
s<R „,  it  readily  follows  that 

lim  4"  f  )'5 i(s)s*rfs  =  -{  •  (2.49) 

R  2*  *  y 21  Jq 

Hence. 


-  nzi  (2  [  1  -  An  Pi  Ci \ (0)  |  -  1 1  —  A',|  pt  C|i(0)  |  | 


(2.50) 


by  inverting  Kq.  (2.23)  and  (2.24)  and  using  Eq.  (2.27,. 
Hotli  <Vj.(.  )  (for  the  closure  rules  considered  in  the  next 
see  turn)  and  )  are  short-ranged  functions  and  Una 
to  zero  .is  .  -  •  .  and  the  asymptotic  lorm  of  /if , (c >  1-, 
therein  1 

Inn  //f.i  1  3A,,  .  12.72 

which,  follows  iron!  Eq .  (2.27)  and  the  tact  that  /if.1' 
also  short  r  eined.  The  asymptotic  lorm  ol  the  w.iii- 
particle  corn  latum  function  lor  the  HLHNC  closure, 
discussed  111  llu  next  section,  coincides  with  that  ;  :  in 
MS  approximation''  and  is  therelore  given  by 


Inn  //;,(  .  E;.  !.’,)  3A',,Z/(2.  1) 

(2.  33' 

3m,  *3/7(2.  D 

a*/  . 2^ 2A  1  i  > >  1  1  ] )  -  *“  A  1  j  /  )]  Ajj)  1 

'2.  54) 

3FSM2.  1)| 

'  T*  i2r^2K;, ,),!([,)  .liS-K,,., ,/{[,)';  ’ 

(2.  5' ' 

where 

(2.  at 

and 

0  //.  j  (3  cosafr2  -  1)  '  ‘  n/  j 

(2. 57’ 

Since  the  electric  field  is  independent  of  the  distance 
from  the  wall.  h2 E2.  S^)  does  not  decrease  to  zero 
as  ■?  —  «<•,  unless  the  magnitude  of  the  electric  field  or 
/7(2,  1)  is  also  zero.  The  angular  average  ol  /i?1  ( v .  E:.  S.’ 
is  however  zero  in  the  limit  <-K: 

lim  /  /1 21  ( v .  e2,  n,)r/n,  =  o 

e- -  J 

(2.58) 

since  the  angular  averages  of  £>(2,  1)  and  A(2,  1) 
and 

are  zero 

lim  h^U)  =  0  . 

(2.  59) 

We  shall  now  consider  the  details  of  these  equations  in 
the  context  of  the  mean  spherical  and  linearized  hyper - 
netted  chain  closure  rules. 

HI.  THE  CLOSURE  RELATIONS 

The  closures  for  the  wall -particle  direct  correlation 
functions  arc  readily  derived  by  taking  the  wall  limit  ol 
the  corresponding  closures  in  the  bulk  fluid. 


J  Chem.  Phys.,  Vol.  73.  No  8.  15  October  1980 


3984 


Eggebrecbt,  Isbister,  and  Rasaiah  Dipoles  at  a  wall 


A.  MSA 

Tilt'  closure  tor  the  direr!  rrolatiun  function  is 


The  .ttiovc  equations  13.  121  'Acre  solved  iterate.  '  r. 
(onq.  ili  r  to  "ini'r.ilr  t lit-  accompany nip  ti r'a n  . 


/XI.  2)  .  Ion  .  (3.1) 

A(  y 

which  is  equivalent  to 

r^U-i  0  .  tor  r  •  /<„,  .  (3.  21 

This  is  unaffected,  tor  K:I.  oil  taking  the  wall  limit 
of  cl,(r! . 


B.  LHNC  approximation 

Starting  with  the  hypei  .•lotted  chain  (UNO  approxima- 


rn„(2, 1) -;ia42.  11  -  tnw(2.  1>  -  l\,6(2.  D  AT  (r  •«„) 

(3.3) 

and  using  the  invariant  expansions  of  ra6( 2,  1)  and 
/ia?( 2.1)  ;cf.  Eqs.  (2.131  and  <2. 14)  |.  one  finds 


at 


M2.  1)  . 


where  \'*s  is  the  spherically  symmetric  part  of  the  po¬ 
tential.  Expanding  the  logarithm  up  to  the  linear  term, 
and  collecting  and  comparing  coefficients  of  1.  M2,  1), 
and  -M2,  1).  we  have 

r •£„<»•)  ■-  h‘at[r)  -  ln^«fl(»  )  ~  '•>,(>  )  AT  ,  (3.5) 


r3  AT  ' 


IV.  RESULTS  AND  DISCUSSION 

Tin  solution  of  E<j.  (2.51).  for  A2,(.  .  £•..  l.:'.  a.  -  ■* 
lamed  as  tiie  confluence  of  three  distinct  computation- 


(li  fiie  bulk  correlation  functions  wt  .  calculator; 
tin  niaiiiu-r  described  by  Pa  ley,  from  Eq.  (2.37.  T  - 
reduced  difx.le  moment  dependence  ul  this  relation  . 
tallied  in  llie  constant  A’n,  through  Eq.  (2.38k  i...c. 
lie  expressed  in  the  MSA  and  ULHNC  approximations  a - 


whore  r{',(«J,>  and  6f,(/fJ,)  are  the  values  o!  ci'.fr)  and 
b{.(i  )  immediately  outside  contact.  The  value  o! 
was  adjusted  until  this  difference  assumed  the  desirid 
reduced  riijxde  moment.  The  iteration  ol  Eq.  (2.3. 
was  performed  using  last  Fourier  transform  tevnniqu.  - 
and  mixed  solutions  to  speed  convergence.  Tin  ou!-. 
ituid  correlation  functions  obtained  showed  c.\c< 
aurc-enie.’i!  with  those  ot  Patey." 


(2'  Tlie  electric  field  which  emerges  Iron,  tin  w.d. 
was  determined  by  Eqs.  (2.9),  (2.30).  and  12.  40'-.  j  h. 
field  angle  a,  of  Fig.  1  was  taken,  in  separate  i-.ilr.il..- 
nons.  as  O'.  45',  and  90".  The  field  angle  is  related 
to  the  wall-dipole  orientation  through  Eq.  (2.  KM.  the 
solution  of  winch  appears  in  Fig.  1.  The  magnitude  ol 
the  reduced  electric  field 


E2 :  ft?,  E|,(3 cos2t>g s  l)1/2B2 


<■£.(»■>  =*i,(r)[l  -ATl. (»■)■' : . 

(3.7) 

Defining 

C(>-)  =  C(>')[i-.Aie(>r1i  - 

(3.8) 

*£»(»')  =rt£s(r)  [1  -AW !  - 

(3.9) 

with  h»fl(r)  and  b'ae(r)  also  defined  by  equations  analogous 
to  Eqs.  (2.18),  (2.23),  and  (2.24),  respectively,  the 
last  two  equations  can  be  written  in  the  form 

clsW  =  ftaeW  »  r>Ra „.  (3.10) 

In  the  RLHNC,  Eq.  (3.5)  is  replaced  by  the  exact  closure 
for  hard  spheres  against  a  wall. 

In  the  wall  limit,  6f,(e)  =bfx(z),  and  making  use  of 
Eq.  (2.27),  the  closure  condition  for  the  wall -particle 
direct  correlation  function  becomes 

C|i(^)  =  [l  +/i|iU)J[l  -2etiUr']  s  *>0  .  (3.11) 

On  substituting  this  in  Eq.  (2.34),  we  have,  for  z>0, 

**„<*)  =  [*!,(*>  "  1 1*  *Si(*>  2^u  P\  {*  [«*(-)  -B*M  I 

'+[£>*<«)-£>*(*)  I  +  J  *|,U)lB*(«)-B*(|2-y|)Uv|  • 

(3. 12) 

When  g|,(<-)  - 1,  the  formal  equation  for  /i|[(z)  in  the 
MSA  approximation  is  recovered. 


was  taken,  as  in  the  earlier  work  of  Isbister  and 
Freasicr,  lo  be  with  the  reduced  dipole  momem 
squared  »*!  fixed  at  either  0.5  or  1.0.  For  purpost 
of  comparison  of  these  parameters  to  a  molecular  sy~- 
ti  m,  the  reduced  dipole  moment  squared  for  HC1  (o.: 

■  1.03  l»  at  275  "K.  assuming  a  diameter  of  3.5  A,  is 
approximately  0.65.  A  reduced  field  Strength  o!  B  3 
for  this  system  corresponds  to  an  electric  field  of  near¬ 
ly  1.9  •  109  V'ni  or  a  surface  charge  density  a  of  1  elec¬ 
tronic  charge' 1000  A2. 

Wt  have  also  carried  out  calculations  at  the  same 
surface  charge  density  (or  electric  field  £2)  for  a  fluid 
at  a  reduced  density  of  0.  7  with  the  reduced  dipole  mo¬ 
ment  w*  -  2.0.  These  numbers  correspond  approxi¬ 
mately  to  those  appropriate  for  liquid  water  (hi,  -  1 .  B5 
D,  /?,,  =  2.  76  A)  at  room  temperature.  The  reduced 
electric  field  £|  [which  contains  m,  and  £j,  in  its  defini¬ 
tion  (4.  2)  j  is  now  only  0.71.  (It  may  be  useful  for  the 
reader  to  bear  in  mind  that  it  is  an  artifact  of  our  defini¬ 
tion  of  E*  that  an  increase  in  the  dipole  moment  w,  re¬ 
sults  either  in  a  reduction  of  £j  when  the  surface  charge 
density  is  held  constant,  or  an  increase  in  the  surface 
charge  density  when  Ej  is  unchanged. ) 

(3)  The  spherically  symmetric  part  of  Eq.  (2.  51). 
i.  c. ,  AJ,(.:),  was  determined  using  the  technique  of 
Watsinan.  Henderson,  and  Lebowitz. 7  The  function 
/i|,(z)  was  first  obtained  from  Eq.  (2.28)  as  the  Porcus- 
Yevick  solution  and  then  corrected  to  produce  an  essen- 


J,  Chem.  Phy$.,  Vof  73.  No  8,  1b  October  1980 


Eggebrecht,  Isbister,  and  Rasaiah:  Dipoles  at  a  wall 


986 


k. 


b  f 


R 


5R, 


K1G.  The  wall- particle  distribution  functions  g  \U>  ii.i 
as  a  function  of  the  distance  z  from  the  wall  for  different  orien¬ 
tation#  U«  of  the  fluid  dipoles  where  £*  =  v3,  n\  f  ‘  -  o. p* 

*0.373,  and  o.  »t«  .  /t.  is  the  radius  of  a  fluid  dipole. - 

IM.HNV  approximation, - MS  approximation. 


-2- 


7  ~ 


b  t1, 


5P, 


FIG,  8.  Wall-particle  dis’  ribution  functions  for  the  system 
depicted  in  Fig.  2  except  ttisi  a,  =  4 5". 


FIG.  J.  Wall-particle  distribution  functions  for  the  systems 
depicted  in  Figs.  2  and  :i  except  that  a,  =  9(1°. 


3R 


FIG.  f>.  The  wall-particle  distribution  function  y  ■  .  E.,  1 1,1 

for  the  system  depicted  in  Fig.  I  except  that  m"  l.o. 


I  20  -  '  '  ’  r~T 

I 


0  PC  V- 


)  hi.  The  expansion  coefficients  h*  U)  and  kl\  i.?>  as  a  Junc¬ 
tion  oj  the  wall-particle  distance  z  when  £?  r  h  :i  and  p?  -  11,573, 
Hie  upper  and  lower  curves,  in  each  case,  are  u>r  ».■?  ] .  0 

and  0.5,  respectively. 
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l'iR, 


I'M,  7 .  Hie  w:  ll  —  i >:■  i-lioK-  ,u*u  tier i.m  (unction  g  .  m,  E  ,  I.M 
as  a  hmclion  ot  liu*  distance  tr  *  *  toe  wait,  lor  ditforenl 
orientations  I.’,  oi  tin1  tb;ii  dipole*.  when  E*  ‘*.~l,  ,)I *  t.ii, 
p;  0.  7,  anil  : > .  0.  [’!»•  rcd'iee.i  parameters  m*  anil  ,»7  art* 

those  appropriate  tor  liuuiil  w.ner  .3  room  temperature  on, 
-t.Sa  l),  8tl -*  .’.76  \i.  llu*  reilueetl  electric  field  £*  0.71 

corresponds  to  a  surface  densiti  ot  1  electronic  charge  loot'  A  . 


tiully  exact  wall -particle  correlation  function  in  com¬ 
parison  to  Monte  Carlo  simulation.  The  reduced  density 
of  hard  spheres,  upon  which  h5,(. r)  depends  solely,  was 
taken  as  p*  =0.  573. 

The  computations  were  combined  through  Eqs.  (2.23), 
(2.  24),  and  (3. 12)  to  yield  the  coefficients  of  the  angular 
functions  of  Eq.  (2.51).  The  integrals  over  bulk  corre¬ 
lation  functions  (2.33)  as  well  as  the  integral  of  Eq. 

(3.  12)  were  performed  by  trapezoidal  approximation. 

The  iterative  solution  to  Eq.  (3.12)  was  rapidly  conver¬ 
gent.  Finally,  the  expansion  coefficients  were  inserted 
into  Eq.  (2.51). 

In  Figs.  2-5  we  present  our  results  for  g21(z.  E2,  ST ,) 
in  the  KLHNC  approximation.  The  expansion  coeffi¬ 
cients  h?,(r)  and  /t2l(r)  are  shown  in  Fig.  6.  A  compari¬ 
son  of  KLHNC  and  MSA  results,  for  a  field  angle  of  zero 
and  reduced  dipole  moment  squared  of  0.5  and  1.0,  are 
shown  in  Figs.  2  and  5,  respectively. 

In  the  case  of  m*2  =  0.  5.  F,'  -  8  3,  and  p*  -  0.  573 
tf21(z,  E2,  H,)  in  the  KLHNC  approximation  assumes 
negative  values  only  very  near  and  at  contact  for  a  di¬ 
pole  orientation  in  opposition  to  the  field.  For  higher 
values  of  m*2-  1.0,  at  the  same  reduced  electric  field 
and  fluid  density,  these  regions  are  more  pronounced 
and  extended  in  range  (Fig.  5).  An  increase  in  w*2 
from  0.5  to  1.0  at  constant  E*  and  p* ,  however,  corre¬ 
sponds  to  a  fourfold  increase  in  the  surface  charge  den¬ 
sity  from  1  electronic  charge  1000  A2  to  one  every 


250  A2,  liolh  the  KLHNC  and  MSA  treatments  result 
in  distributions  showing  exhanced  adsorption  lor  favor¬ 
able  di|x>lc  orientations,  but  flu*  repulsive  interaction.- 
of  the  electric  field  with  unfavorably  aligned  dipoles  1- 
niorc  clearly  visible  in  the  KLHNC  approximation.  Tins 
theory,  like  the  MSA.  is  essentially  linear  in  character, 
and  cannot  prevent  the  distribution  functions  from  be¬ 
coming  negative  when  the  dipoles  are  aligned  against  ti.< 
field,  it  ai  the  same  time  the  theory  predicts  a  large  en¬ 
hancement  ot  the  adsorption  of  dipoles  aligned  with  the 
field.  When,  for  example,  tile  electric*  field  is  perpen¬ 
dicular  lo  the  wall,  both  theories  predict  that  the  den¬ 
sity  proliles  (see  Figs.  2  and  5)  are  symmetrical  about 
the  profile  tor  dipoles  perpendicular  to  the  Held  f:  . 

-r.  2  or  3:  21.  A  large  enhancement  of  dipoles  align-  i 
with  the  field  would  lead  to  an  equally  large  depletion  ■* 
dipoles  oriented  against  the  field,  which  may  require 
negative  wall -particle  distribution  functions.  Eurthi  : 
improvements  beyond  the  RLHNC  approximation  so.il>; 
then  be  necessary.  Figure  7  shows,  however,  that  It,* 
KLHNC  theory  provides  plausible  density  proliles  i  ven 
when  the  square  of  the  dipole  moment  n.*‘  is  increased 
to  4.0  hut  the  surface  charge  density  is  maintained  at  1 
electronic  charge' 1000  A2. 
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